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ABSTRACT 

We explore the possibility that large-scale convection be inhibited over some regions of giant planet 
interiors, as a consequence of a gradient of composition inherited either from their formation history 
or from particular events like giant impacts or core erosion during their evolution. Under appropriate 
circumstances, the redistribution of the gradient of molecular weight can lead to double diffusive 
layered or overstable convection. This leads to much less efficient heat transport and compositional 
mixing than large-scale adiabatic convection. We show that this process can explain the abnormally 
large radius of the transit planet HD209458b and similar objects, and may be at play in some giant 
planets, with short-period planets offering the most favorable conditions. Observational signatures of 
this transport mechanism are a large radius and a reduced heat flux output compared with uniformly 
mixed objects. If our suggestion is correct, it bears major consequences on our understanding of giant 
planet formation, structure and evolution, including possibly our own jovian planets. 
Subject headings: convection — hydrodynamics — planets and satellites: general — stars: planetary 
systems 
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1. INTRODUCTION 

Over 200 extrasolar giant planets have now been dis- 
covered by radial velocity surveys. Fourteen of these 
planets have been observed transiting their parent star, 
allowing an accurate determination of their radius and 
mean density. For half of these objects, notably the first 
ever discovered transit HD209458b, the predicted theo- 
retical radius lies sev eral cr's (up to ~ 20 %) below the 
observed mean value (jBaraffe et al. 20051 ). The various 
scenarios proposed so far to solve this discrepancy either 
have b een rejected on obser vational or theoretical argu- 
ments (jLevrard et al. 20 07) or lack an identified robust 
mechanism to conver t surface kinetic energy in to ther- 
mal energy at depth (|Showman fc Guillot 20021) . A fair 
conclusion of these studies is that an important physical 
mechanism is probably missing in our present description 
of at least some short period planets, for which we have 
a radius determination, but possibly of all extrasolar or 
even solar giant planets. 

According to the conventional core-a ccretion model for 
planet formation (Pol lack et al. 19961 ). planets are be- 
lieved to have a substantial enrichment in heavy elements 
compared with their pare nt star, with a total of ~ AOMq 
for a Jupiter-mass object (|Alibert et al. 20 05) . Observa- 
tional constraints of Jupiter and Saturn show that these 
planets do have a significantly enhanced Z-abundance 
compared with the Sun, w ith a global mean mass frac- 
tion Z=M Z /M p « 10-20% dSaumon fc Guillot 20041) . In 
all these calculations, the big planetesimals are gener- 
ally supposed to drown to the core during the early 
phase of solid accretion, while the smaller ones are dis- 
tributed uniformly throughout the envelope, leading to 
a uniform heavy element abundance. This is a rather 
simplistic description of the planet internal structure, 
which implies (i) well defined interfaces between the 
central core and the (highly diffusive) H/He rich enve- 
lope and (ii) very efficient large-scale thermal convection 



throughout the entire gaseous envelope. The observed 
atmospheric abundances of Jupiter and Saturn, however, 
seem to require the redistribution of a subsequent frac- 
tion of heavy eleme nts in the interior of these planets 
(|Guillot et al. 2004). In this Letter, we explore the con- 
sequences of the presence of an initial compositional gra- 
dient in the envelope, as a result of either early plan- 
etesimal accretion or subsequent core erosion, and of the 
resulting less efficient heat transport and compositional 
mixing on the fate of gaseous planets. We show that lay- 
ered convection, if it occurs as a result of this composi- 
tional gradient, might be the lacking physical mechanism 
to explain the transiting planet abnormally large radii. 

2. LAYERED SEMICONVECTION 
The presence of a positive compositional gradient, i.e. 



a gradient of mean molecular weight V 
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tends to stabilize the fluid against convective instability 
according to the Ledoux stability condition: 



XT 



(1) 



where Vt and ~S7 a d denote the usual temperature and 
adiabatic gradients, respectively, and Xm = ( 3in^ )p- T ' 

XT = ( g l £ lT ' ) p<v ^ n mos t °f the giant planet interiors, 
superadiabaticity is extremely small, with Vt — V a d ^ 
10~ 8 , so a small molecular weight gradient over a typ- 
ical mixing length size region can affect significantly 
and even damp out convection. In convective systems 
where buoyancy effects of (destabilizing) heat and (sta- 
bilizing) composition are opposed, the process leads gen- 
erally to quasi-static uniformly mixed convective layers 
separated by small diffusive interfaces with steep gradi- 
ents, V M 3> 1. Such stable layered convection is indeed 
observed in some areas of the Earth's oceans, due to 
the presence of the stabilizing salt gradient (thermoha- 
line convection) leading to a stratified ste p-like tempera - 



ime convection ) leading to a stratmed ste p-iiKB tempera - 

1 and Max Planck Institut for Astrophysics, Garching, Germany ture profile with stable boundary layers (jSchmitt 19941) . 
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Lab oratory experim ents have also confirmed this layer- 
ing (jFernando 19891 ). Once formed, the stratification is 
stable provided the compositional gradient remains large 
enough at the interfaces to satisfy Eq.([T]). Would such 
a stratification occur in planetary interiors, the layered 
part of the interior can be considered as a semiconvec- 
tion zone with a reduced efficiency to transport the inter- 
nal heat and composition flux compared with large-scale 
convection. One may argue that the conditions in plane- 
tary interiors differ from the ones in the oceans or in the 
experiments. Characteristic thermal diffusivity in H/Hc 
planetary interiors, dominated by electronic transport in 
the central, ionized parts, and by molecular motions in 
the outer envelope, lies in the range kt ~ 10 -1 -10~ 2 
cm 2 s _1 , while the kinematic viscosity is v w 1CP 3 -10 -2 
cm 2 s _1 ()Stevenson fe Salpeter 1977al ). The characteris- 
tic Prandtl number thus ranges from Pr = v/kt w 10 -2 
to 1. These values do not differ by large factors from the 
ones characteristic in the oceans or in laboratory exper- 
iments, Pr w 1-10, in contrast to the ones characteristic 
of stellar conditions (Pr < 10 -6 ). Therefore, a layering 
process in giant planet interiors can not be excluded and 
it is worth exploring the consequences on the planet evo- 
lution. Such layered convection may occur near the dis- 
continuity in composition at the boundary of the central 
rocky-icy core or in chemically inhomogeneous regions in 
the interior, reminiscent of the early planetesimal accre- 
tion episodes. 

There is presently no widely accepted descrip- 
tion of semic onvection. Water-salt experiments 
((Fernando 19891 ) show that a series of quasi-static con- 
vective layers separated by diffusive interfaces develop 
when a balance is reached between the variation of po- 
tential energy (i.e. of buoyancy) due to mixing at the 
interface and the kinetic energy of the eddies available 
at the interface. This translates into a critical Richard- 
son number Ri—lgAp/p(v 2 ) of order 1 to 10, where 
Ap ~ p(AT/T) is the density contrast between the dif- 
fusive and convective layers, g is the gravity and p(v 2 ) is 
the kinetic energy of the convective flow, of characteristic 
average length scale I and rms velocity \J (v 2 ). G uided by 
experimental results ([Linden fe S hirtcliff e 1978T ) and en- 
ergetic arguments, Stevenson (1979), in a wave descrip- 
tion of semiconvection, showed that, should layers form 
as a result of small-scale wave breaking whereupon the 
compositional gradient is redistributed, they would be 
stable if 

( KT + ^)(V T - V ad ) < (I> + i/)V M , i.e. Pr>T l ' 2 {2) 

where the inverse Lewis number T—Le^ 1 —D/ kt is the 
ratio of the solute microscopic diffusivity to the thermal 
diffusivity. Spruit's (1992) stability condition is less re- 
strictive, since layered formation is supposed to always 
occur, and is given essentially by Eq.(l). Under jovian 
planet c onditions, typical values are D w 10~ 3 -10~ 4 
cm 2 s _1 (jStevenson fc Salpeter 1977al) . then r - 10~ 2 , 
so that, according to this criterion, diffusive layers could 
be at least marginally stable in giant planet interiors. It 
is worth noting that the molecular to thermal diffusiv- 
ity ratio is the same for a H/He mixture under jovian 
interior conditions as for salty water, r k, 10~ 2 , so that 
the extent of the solute versus thermal layer is about 
the same. Since, according to experiments and condition 



(2), this ratio is the relevant criterion for stability of the 
layers, this adds some support to the planetary case. 

In a layered convection stratification, heat is carried 
away from the interfaces by descending and ascending 
plumes in the overturning regions while transport across 
the interface occurs by diffusion. Because of the bound- 
ary layers, only a part of the fluid transports heat ef- 
ficiently. On average, one has "convective- like" motions 
having a much shorter length scale than for ordinary con- 
vection. The thermal thickness of the diffusive layers, St, 
is determined by a balance between the thickening due 
to diffusion and the entrainment due to convective mo- 
tions so the convective time, ~ l/v, in the mixed layer of 
size I must be comparable to the therm al diffusion time , 
~ S t /kt, across the boundary layer (jFernando 19891) . 
This yields 5 T w (k t l/v) 1 / 2 and S x » (Dl/v) 1 / 2 for the 
thickness of the heat and compositional interfacial layers, 
respectively. The number N of layers is of course very 
uncertain. It can be crudely estimated as follows. The 
heat flux / transported by convection in each mixed layer 
is the mass flux carried by the plumes fed from the diffu- 
sive layers, and thus of width comparable to these layer 
thickness, times the energy variation across the convec- 
tive layer, TSS « T(l VS) : 

f = ]( P vS T ) x (TSS) = pcpKT-^^W - V ad ) (3) 

For a semiconvective region extending over a planet-size 
region, the total number of layers is thus given by: 

N ^^i = L J^ X *L (4 ) 

/ pc p K T jf^(S7 - V ad ) I 

Using characteristic numbers for jovian planet conditions 
(with / ~ Hp), one gets N ~ 10 5 -10 6 as a rough estimate. 
If convection is inhibited, however, the smaller heat flux 
and larger superadiabaticity require less layers. 

3. EFFECT ON THE EVOLUTION 

We have conducted calculations following the evolu- 
tion of a template Jupiter-mass planet, representative 
of HD209458b and similar short-period planets, with 
a global metal content Mz — 40 M®, including a 6 
M® core, i.e. M z /M p = 13% (Z ~ 6Z Q ), in agree- 
ment with previously mentioned planet formation mod- 
els and Jupiter and Saturn's observational constraints. 
The amount of heavy elements is distributed initially 
throughout the planet following a gradient, distributed 
within a certain number of boundary layers N where con- 
dition (1) is fulfilled. The layers are located within the 
inner ~ 30% by mass (60% in radius) of the planet, where 
H and He are fully ionized, to ensure high enough ther- 
mal conductivity. The present calculations have been 
done with N=50 and 7V=100; the width of each bound- 
ary layer corresponds to St ~ 10 3 cm ~ 10~ 6 iJp. A 
larger number of layers would be computationally too 
difficult to resolve correctly. These boundary layers are 
separated by larger convective, mixed layers, with a uni- 
form composition (V/i = 0), where the usual mixing- 
length formalism applies. The sizes of the boundary and 
mixed layers (St,Sx,1) obey the aforementioned rela- 
tionships. The heat flux Ft and solute flux Fz in the 
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boundary layers are calculated with the appropriate dif- 
fusion equations : Ft = pc p ktVT, where kt = ^§p^T ( K c 1 + 
K r ^) and k C: K r denote the conductive (Potckhin 1999) 
and radiative (|Rogers et al. 1996t iFerguson et al. 2005ft 
mean opacities, respectively, and Fz~pD\7Z, where 
VZ=-^Hp l V lli with x = (mtlz )> is the mean concen- 
tration gradient across the layer. Conduction remains ef- 
ficient enough in the thin boundary layers to fulfill condi- 
tion (1). Because diffusion limits the heat transport, the 
internal heat flow of the planet is significantly reduced 
compared with that of a fully convective object. The 
signatures of double-diffusive convection in a planetary 
interior are thus a reduced heat output and a larger ra- 
dius compared with an object where heat is transported 
efficiently by large-scale convection. This is illustrated 
in Fig. 1, which compares the evolution of the radius 
and thermal intrinsic luminosity of the planet in both 
cases. The excellent agreement with the otherwise unex- 
plained observed radii of HD209458b and similar irradi- 
ated planets suggests that diffusive convection might be 
taking place in the interior of at least certain giant plan- 
ets. As seen, the expected luminosity at young ages is 
more than one order of magnitude fainter than that of a 
fully convective planet evolving from a comparable initial 
state. The observational confirmation of the present sce- 
nario would be either the determination of an exoplanet 
temperature or luminosity at young ages 2 or the observa- 
tion of an inflated radius for a transiting planet at large 
enough orbital distance, a>0.1 AU for a solar-type par- 
ent star, for stellar irradiation not to affect the planet's 
internal structure. Figure 1 illustrates also the depen- 
dence of the evolution upon the number of layers. Less 
boundary layers implies larger convective layers and thus 
more efficient heat transport, as illustrated by the more 
rapidly decreasing radius in the 50-layer calculations. 

A key question is to know if diffusive interfaces can 
persist on time scales comparable to the characteris- 
tic time for the evolution of the planet. According to 
the aforementioned critical Richardson num ber criterion, 
supported by experiments ([Fernando 19891 1 . a quantita- 
tive argument is that if the average kinetic energy in the 
convective layers is smaller than (a fraction of) the po- 
tential energy wall of the interface, convection cannot 
penetrate deeply into this latter and significant entrap- 
ment across the interface cannot occur. This implies that 
the molecular diffusion time scale be long enough. This 
latter can be estimated for the entire stack of layers, dis- 
tributed over a region of size L in the planet (presently 
L ~ 10 9 cm). The flux of element across an interface 
is Fz ~ pDj^, where SZ « j- AZ is the jump in the 
element mass fraction at each interface while AZ is the 
total variation over the entire semiconvective region. The 
time scale to redistribute the entire gradient over the en- 
tire region is then t w w ^ w 10 Gyr. This 
admitedly crude estimate shows that the stable diffusive 
convection configuration might last long enough to af- 

2 For short-period, irradiated planets, however, the intrinsic flux 
of the planet, crT^.p is smaller than the absorbed and reflected 

contributions of the incident stellar flux, ~ (i?*/a) 2 jF*. For long- 
period planets or for planed telescopes like the LBT or the JWST, 
dedicated to infrared planet searches, however, the planet intrinsic 
luminosity can be determined. 
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Fig. 1. — Evolution of the radius and the intrinsic luminosity 
of a Jupiter-mass (1 Mj up ) planet orbiting at 0.05 AU of a Sun- 
like star, from the same initial conditions. All calculations include 
the effect of the stellar irradiation on the planet structure and 
evolution. (1) Effect of a core and metal enrichment: dotted line: 
adiabatic interior, no core, Z = Zq; long-dashed line : adiabatic 
interior, central dunite core M c = 6M^, metal enrichment Z env ~ 
6 X Zq in the envelope. Effect of layered convection (same M c and 
Z e nv)- solid line: layered convection for ./V=100 layers; dot-dash 
line (upper panel): 7V=50 layers. The observed values of the seven 
planets with abnormally large radii, HD209458b (triangle), WASP- 
lb, XO-lb, OGLE-TR-56b, Tres-2, HAT-P-lb and HD189733b, 
with masses ranging from 0.5 to 1.3 Mj up , are displayed with their 
most recent l-cr error bar determinations. 



feet substantially the evolution. With the typical value 
D = 10" 3 cm 2 s" 1 , about 10% of the initial gradient AZ 
has been transported by diffusion over a Gyr, as con- 
firmed by our numerical calculations. In principle, the 
compositional gradient thus remains large enough during 
the evolution for the Ledoux criterion to remain valid in 
a majority of layers. In other words, the temperature 
jump at interfaces is too small to offset the molecular 
weight stabilization of interfaces (^?- < — ). The com- 
position and temperature profiles in our calculations at 
5 Gyr are portrayed in Fig. 2, with SZ/AZ « 1% and 
AT w 10 3 K at each diffusive interface. Note that, if lay- 
ers form in sequence through turbulent entrainment or 
from sporadically breaking internal waves generated by 
oscillatory instabilities, interfaces may be dynamically 
renewed with time, if some compositional gradient or 
stirring effects remain present. Such a process occurs 
in laboratory systems and oceans. 

Different reasons can be advocated for the cause of 
the initial compositional gradient. This latter can be 
inherited from the formation process. Large incom- 
ing planetesimals could disseminate part of their con- 
stituents, iron, silicates, ices, by ablation and break- 
up as they penetrate the building gaseous envelope 
(|Iaroslavitz fc Podol ak 2007). Note also that accretion 
will not proceed homogeneously as capture mechanisms 
differ for the gas (H, He), ice (essentially C, O, N) and 
rock (silicates and iron) components. This will increase 
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Fig. 2. — Internal heavy element and thermal profiles at 5 Gyr 
for the case of 100 convective + diffusive layers distributed in the 
inner 30% by mass of the planet. The global heavy element mass 
fraction of the planet is Z = Mz/M p = 0.13, including the central 
core (Mcore = 2%M P ). Note the step-like T and Z-distributions, 
as portrayed in the inner subsets. 



substantially the compositional gradients before the core 
is reached. As mentioned earlier, even modest gradients 
can easily offset superadiabatic excess over planet-size re- 
gions, preventing large-scale convective motions. A com- 
positional gradient might also result from disruption and 
redistribution of the core due to a giant impact or erosion 
at the core-envelope interface because of metallic hydro- 
gen high diffusivity . leading to a co re diluted into a frac- 
tion of the planet (|Stevenson 19821 : iGuillot et al. 2004t ). 
The redistribution of these elements might be partially 
inhibited by diffusive processes, forming diffusive inter- 
faces because of the opposite buoyancy effects of heat 
and composition. Furthermore, when accreting the en- 
velope, part of the outermo st regions of th e protoplanct 
might be nearly isothermal (jMizuno 1980f ), which favors 
the stability of a compositional gradient. Interestingly, 
when distributing the layers in the outer 10% by mass 
(~ 15% by radius) of the planet, where radiative thermal 
diffusivity starts to dominate over conductive diffusiv- 
ity, we get an effect similar to the one portrayed in Fig. 
1. The situation for the formation of diffusive convec- 
tion is particularly favorable for short-period exoplanets 
for several reasons. First of all, a substantial fraction of 
the gaseous envelope has been eliminated by evaporation 
(|Baraffe et al. 20061 ) . leading to a larger metal fraction. 
Second of all, for short-period exoplanets, the numerous 
collisions tend to eject the gas, leading to a larger enrich- 
ment in planetesimals than for the other planets. Third 
of all, the higher internal temperatures for short-period, 
irradiated planets, favor (i) ionization of the various el- 
ements and thus the thermal conductivity, (ii) solubility 
of the core material into the envelope. At last, because 
of the stellar irradiation, the outer layers of short-period 
planets are isothermal and not adiabatic. 

4. CONCLUSION AND PERSPECTIVE 

The aim of the present Letter is to suggest an alter- 
native, possibly important energy transport mechanism 
in giant planet interiors and to explore its effects on the 



evolution. These calculations provide a consistent de- 
scription of the evolution of giant planets, with a metal 
enrichment in agreement with observational constraints, 
in case heat is transported by layered convection. Assum- 
ing an initial compositional stratification within a cer- 
tain number of double-diffusive interfaces, diffusive and 
convective transport in the respective layers are calcu- 
lated consistently during the evolution. Only the out- 
ermost 40% in radius (70% in mass) of the planet can 
convect freely. These calculations, however, cannot be 
expected to give an accurate description of the onset 
and stability of layered convection. There is presently 
no accurate treatment of this mechanism under condi- 
tions characteristic of giant planets. The only attempt 
to study the onset of double-diffus ive layer forma tion 
at low Prandtl numbers (Pr < 1) (jMerryfield 19951 ) re- 
mains inconclusive. Indeed, insufficient numerical res- 
olution and artificially enhanced viscous and molecu- 
lar diffusivities in the simulations might suppress small- 
scale motions/instabilities and the formation of a statis- 
tically steady state of intermitent diffusive layers. Simu- 
lations of vertical salinity in water, Pr=7, on the other 
hand, well reproduce the experiments and confirm the 
formation of quasi-static convective layers separated by 
diffusive interfaces above a c ritical Richardson number 
(jMolemaker fc Diikstra 1987T ). 

Even though the present calculations rely on some un- 
certain ground, the agreement with the puzzling and oth- 
erwise unexplained observed radius of HD209458b and 
other abnormally large exoplanets leads to the conclusion 
that layered convection might be taking place in at least 
some planets and could explain their particular proper- 
ties. This transport mechanism yields a much reduced 
heat escaping rate compared with a homogeneous adi- 
abatic structure. Even if a stable layered configuration 
does not occur, overstable modes of convection (fullfilling 
Eq.(l) but not Eq.(2)), due to the presence of opposite 
diffusive processes (composition and heat) of different ef- 
ficiencie s, can lead to the growth of sm a ll-scale fluid osci l- 
lations (jStevenson fc Salpeter 19775 IStevenson 1979h . 
Overstability, however, is more similar to an enhanced 
diffusion process than to a convective mixing process, 
with a much sma ller energy transport efficiency , as shown 
by experiments (jStevenson fc Salpeter 1977bf i. The on- 
set or persistence of layered or overstable (oscillatory) 
convection might require optimal conditions, inherited 
from particularly favorable formation or evolution his- 
tories (e.g. late accretion of large unmixed planetesi- 
mals or giant impacts stirring up completely the planet 
interior). The present paper suggests that, under such 
appropriate conditions, the heat transport mechanism in 
giant planet interiors can be severely affected, decreasing 
the efficiency of or even inhibiting large-scale convection. 
This should motivate 3D investigations of convection in 
the presence of a stabilizing compositional gradient un- 
der conditions suitable to giant planets and the search 
for transits at larger orbital distances. 

The authors are indebted to the MPA for the warm 
hospitality. We are grateful to Henk Spruit for valuable 
discussions and to Friedrich Kupka for pointing out use- 
ful references. 
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Fig. 3. — Evolution of the radius and the intrinsic luminosity of a Jupiter-mass (1 Mj up ) planet orbiting at 0.05 AU of a Sun-like star, 
from the same initial conditions. All calculations include the effect of the stellar irradiation on the planet structure and evolution. (1) 
Effect of a core and metal enrichment: dotted line: adiabatic interior, no core, Z = Zq; long-dashed line : adiabatic interior, central dunitc 
core M c = 6M^, metal enrichment Z env ~ 6 X Zq in the envelope. Effect of layered convection (same M c and Z env ): solid line: layered 
convection for 7V=100 layers; dot-dash line (upper panel): 7V=50 layers. The observed values of the seven planets with abnormally large 
radii, HD209458b (triangle), WASP-lb, XO-lb, OGLE-TR-56b, Tres-2, HAT-P-lb and HD189733b, with masses ranging from 0.5 to 1.3 
Mj up , are displayed with their most recent 1-a error bar determinations. 
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Fig. 4. — Internal heavy element and thermal profiles at 5 Gyr for the case of 100 convective + diffusive layers distributed in the inner 
30% by mass of the planet. The global heavy element mass fraction of the planet is Z = Mz/M p = 0.13, including the central core 
(M core = 2%M P ). Note the step-like T and Z-distributions, as portrayed in the inner subsets. 



